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Interconnected network of transmission and distribution lines lead to losses 
in the system and weakening the voltage stability in the system. Installing 
Static VAR Compensator (SVC) in power system has known to improve the 
system by minimizing the total loss and improve the voltage profile of the 
system. This paper presents the application of Multistage Artificial Immune 
System (MAIS) technique to determine optimal size of SVC. The 
performance of this technique is tested on the IEEE 14-Bus Reliability Test 
System (RTS). The optimization results show that the proposed Multistage 
Artificial Immune System (MAIS) technique gives better solution of SVC 
compensator planning problem compared to single stage Artificial Immune 


System (AIS) in terms of lower total system loss and improved minimum 
voltage magnitude. 
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1. INTRODUCTION 

Electrical power systems are complex systems consisting of large number of generating units and 
interconnected network of transmission and distribution lines. Furthermore, the power demand is increasing 
at an alarming rate day by day, as a result of which the utilities are being compelled to operate their 
generating units and transmission lines at their maximum capacity [1]. In order to improve the power system 
operation, Flexible AC Transmission Systems (FACTS) devices are used in recent years. Reactive power 
changes produced by load variations and line switching can cause adverse effects on system voltage stability 
and the interconnected system security [2]. 

There are several types of FACTS devices such as Static Synchronous Compensator (STATCOM), 
Static VAR Compensator (SVC), Unified Power Flow Controller (UPFC), Convertible Series Compensator 
(CSC), Inter-phase Power Flow Controller (IPFC), Static Synchronous Series Controller (SSSC), Thyristor 
Controlled Series Compensator (TCSC), Thyristor Controlled Phase Shifting Transformer (TCPST) and 
Super Conducting Magnetic Energy Storage (SMES) [3]. Among the types of FACTS devices, SVC is the 
prominently used by the Grid System Operators (GSOs). One of the reasons for this is SVC is cheaper than 
other FACTS devices especially STATCOM and TCSC. However, SVC is still very expensive in price thus 
its optimal location and size in the network must be optimally ascertained [3]. 

For the past decade, researchers and electrical engineers have introduced many optimization 
techniques for the SVC planning. These techniques include mathematical methods like linear programming, 
quadratic programming and heuristic methods like Evolutionary Programming (EP), Particle Swarm 
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Optimization (PSO) [4] and Artificial Immune System (AIS). AIS is a computational paradigm that belong to 
the computational intelligence family inspired by the biological immune system [5]. The AIS uses an idea 
gleaned from immunology in order to develop systems capable of performing different tasks in various areas 
of research [6]. AIS is a diverse area of research that attempts to bridge immunology and engineering 
developed through the application of techniques such as mathematical and computational modeling of 
immunology. AIS has proven capable of solving many power system problems such as economic dispatch 
problem, distributed generation placement problem and reactive power planning. 

This paper proposes to use Multistage Artificial Immune System (MAIS) technique to solve SVC 
planning problem. MAIS is a new technique that developed from the modification of the original AIS 
technique. Different from the original AIS, MAIS has multiple stages of the same processes of cloning, 
mutation and selection. The purpose of this modification is to make the possible solution broader than the 
original AIS technique. Subsequently, this will increase the accuracy of the results produced as the number of 
comparison increased. The implementation of MAIS to solve SVC planning problem is explained in details 
in the next section. 

It is found that the proposed MAIS is superior than the original AIS to solve SVC planning problem. 
MAIS gives better solution in terms of lower total system loss and improved minimum voltage magnitude 
compared to original AIS. 


2. PROBLEM FORMULATION 

Prior to the employment of MAIS, the optimal location of the SVC units was found using voltage 
deviation index. Voltage deviation index (VDI) is an indicator used to determine how much the bus voltage 
magnitude has deviated from its reference value. A power system is considered to have good voltage profile 
if it has low voltage deviation index value. The formula of voltage deviation index is as shown in equation 


(1). 
St rea A 
AV, = > [== ) (1) 
ii Vref, 
where V,e; is the reference voltage at i™ bus [8] 


The search for the optimal location and size is based on two objective functions, which are to 
minimize total system loss and to improve minimum voltage magnitude. This study considers single 
objective optimization, which means only one objective function is optimized at a time. 


a. Objective Function 1: To minimize total system loss 
The total system loss can be calculated using the following equation: 


br b ob 
Pross = > Rul? = > [V2 + VP = 2ViV; cos; — 5))] Yr COS Vij 2) 

(=0 i=1 j=1 
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where br and b are the number of lines and buses, respectively; R; is the resistance of line 1; His 
the current through line 1"; V; and 5; are the voltage magnitude and angle at node i™, respectively; Y;; 
and ;; are the magnitude and angle of the line admittance between bus i‘ and bus j‘", respectively [8]. 


b. Objective Function 2: To improve minimum voltage magnitude 
Vinin(after SVC) 2 Vinin (before SVC) (3) 


For this second objective function, it is ensured that the minimum voltage magnitude of the system 
is improved after the installation of SVCs in the power system as shown in equation (3). 


c. Operational Constraints 
The two objective functions are subjected to the following operational constraints: 
2.3.1 Real and Reactive Power Balance Contraints 
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The balance of active and reactive powers must be satisfied at each bus. Real and reactive power 
balance with respect to a bus can be formulated as: 


Poi — Pui = Vi > [Vil@%; cos (8: — 5)) + BY, sin(6; — §)]] (4) 
Fi 

Qi — Qui = ay [viler, sin(5; — 6;) — Gj sin(6; — 5))]| (3) 
j=1 


whereP,; and Qg; are the generated active and reactive powers, respectively; P,; and Q,; are the 
load active and reactive powers at nodei, respectively. The conductance, G’;, and susceptance, B’ jx 
represent the real and imaginary components of element Y';; of the [Y’,,] matrix, respectively, obtained by 
modifying the initial nodal admittances matrix when introducing the SVC [10]. 


2.3.2 Voltage Stability Constraint 
The voltage magnitude at all busses are ensured to be within the limits of +5% of its nominal value. 





Vimin s Vinom s Vimax (6) 


3. IMPLEMENTATION OF MAIS TO FIND THE SIZE OF SVCS 

The location of SVC units are determined by using voltage deviation index before the MAIS is 
employed to find the optimal size of the SVC units. For the IEEE 14-Bus Reliability Test System (RTS) as 
shown in Figure 1, it was decided to install only three SVC units based on its size. The locations of the SVC 
units are at buses that have highest VDI values. 








Figure 1. 14-Bus Reliability Test System [9] 


Table 1. Voltage Deviation Index at Base Condition 
Voltage deviation 








Bus Index at each bus 
1 0.0036 
2 0.0020 
3 0.0001 
4 0.0002 
5 0.0006 
6 0.0049 
7 0.0031 
8 0.0081 
9 0.0097 
10 0.0097 
11 0.0066 
12 0.0034 
13 0.0032 
14 0.0039 
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The process of MAIS to find the optimal size of the three SVC units is shown in Figure 2. Firstly, 
three initial values were generated randomly. These three values represent three control variables of SVC 
units to be installed. The series of random number, x; is generated by uniform distribution number [10]: 


x,=rand(a,b)*off_set +min_set (7) 
where; 

1 =1,2and3 

a = number of row 

b = number of column 


off set = maximum range number that can be achieved 
min set = minimum range number 


The random number generated using the uniform distribution number will be assigned as the size of 
the SVC units. Since it was decided to have three SVCs installed, thus three variables namely x1, x2 and x3 
are generated to represent the size of SVC at bus 8, 9 and 10 respectively. The range of random number of 
the SVC units’ size is between -100 MVAR to 100 MVAR. The SVC units can either injecting or absorbing 
the reactive power to or from the buses respectively. 

Equations (3), (4), (5) and (6) are used to ensure the random numbers generated abide the problem 
constraints. Ttwenty individuals of control variables are generated based on the calculated fitness, which can 
be total system loss or minimum voltage magnitude. These individuals are stored in the population pool. 
Subsequently, the twenty individuals are cloned by ten. 

Unlike the single stage AIS, MAIS clones the population with multiple stages. Before the cloning 
process, the initial population are decomposed into four sub-populations, which contains five individuals. 
After that, the sub-populations will undergo cloning, mutation and selection processes. The MATLAB 
programming code for cloning process is as shown in equation (8). 


Clone_k = repmat (A, [i,j ]) (8) 


where; 
k = number of stage 
A = subject to clone 
i = clone the row 
J = clone the column 


For the mutation process, Gaussian Mutation technique is used to produce the off spring. The 
Gaussian Mutation equation is as shown in equation (9). 


Xitmj = ij +N (0.06 ae imin) (4) (9) 
max 
where; 
Xi+mj  —=mutated parent (offspring) 
Xij = parents 
N = Gaussian random variable with mean p and variance y” 
B = mutation scale, 0 <B< 1 
Xjmax = Maximum random number for every variable 
Xjmin = Minimum random number for every variable 
fi = fitness for the i random number 
fmax | =Maximum fitness 


Multistage artificial immune system for static VAR compensator planning (S. C. Mohd Nasir) 


350 qo ISSN: 2502-4752 
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Figure 2. Flow Chart for MAIS Technique 





At the end of each stage, five best individuals are produced made up of twenty individuals. The best 
twenty individuals are then undergo convergence test to find the global optima. 

Convergence criterion is specified when the value of difference between the 20" and the 1“ 
individual is less than 0.00001 for each fitness. Mathematically, the convergence criterion of fitness of total 
system loss and minimum bus voltage are as shown in equation (10) and equation (11) respectively. 


10 
Pross,max a Pioss,min < 0.0001 ( ) 


AVemax _ AVs min < 0.0001 (1 1) 


4. RESULTS AND DISCUSSION 

As mentioned in the previous section, there are two objective functions tested for the MAIS 
optimization. The objective functions are to minimize total system loss and to improve minimum bus voltage. 
Therefore, the results of MAIS to determine the size of the installed SVC units are divided into two parts 
based on its objective functions. For these two objective functions, the locations of the SVC units are at buses 
8, 9 and 10. These locations are determined from the VDI values of all buses where the three SVC units are 
placed at the buses that have highest VDI values. 


a. Minimizing Total System Loss as the Objective Function 
The results of total system loss and minimum voltage magnitude produced before and after the 


optimization for the first objective function are tabulated in Table 2. 


Table 2. Values of Total System Loss and Minimum Voltage Magnitude for Minimizing the Total System 








Loss 
Minimum 
Technique eas Cai: Voltage 
Magnitude (p.u) 
Pre-optimization (without SVC) 13.56 1.0100 
Post-optimization (using AIS 13.48 1.0100 
technique) : i 
Post-optimization (using MAIS 13.46 1.0093 
technique) 





Based on the Table 2, MAIS produced the lowest total system loss compared to pre-optimization 
and AIS results, which is 13.46 MW. However, for this objective function, the minimum voltage magnitude 
produced by MAIS is slidely lower than pre-optimization and AIS results. This value is still acceptable as the 
most minimum value allowed is 0.95 p.u. It is found that SVC units installation can help to minimize the total 
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system loss and by using MAIS technique to find the optimal size of the SVC units will lead to better results 
compared to AIS technique. 


Table 3. Locations and Sizes of SVC Units Based on Objective Function of Minimizing the Total System 











Loss 
Technique SVC location 
Size of SVC Bus 8 Bus 9 Bus 10 
(MVAR) AIS 34.35 35.25 2.79 
MAIS -25.83 32.28 15.29 





It can be seen from Table 3. that all the three SVC units supplied reactive power to the system. 
However, for MAIS, only SVCs at bus 9 and bus 10 supplied reactive power and SVC at bus 8 absorbed 
reactive power of 25.83 MVAR from the system. Different technique leads to different size of SVC units. 


Table 4. Values of Total System Loss and Minimum Voltage Magnitude for Improving the Minimum 
Voltage Magnitude 


Technique Total System Minimum Voltage 
Loss (MW) Magnitude (p.u) 
Pre-optimization (without SVC) 13.56 1.0100 
Post-optimization (using AIS 13.49 1.0150 
technique) : . 
Post-optimization (using MAIS 13.49 1.0200 
technique) 





Based on the results in Table 4. the minimum voltage magnitude of the system improved from 1.010 
p.u before optimization to 1.015 p.u and 1.020 p.u after it is optimized using AIS and MAIS respectively. 
The total system loss also can be seen reduced from the pre-optimization. From Table 5. it can be seen that 
both AIS and MAIS have SVC units that absorbing the reactive power. For AIS, SVC unit at bus 8 absorbed 
the reactive power while for MAIS at bus 10. The other two SVC units of each technique are supplying 
reactive power to the system. 


Table 5. Locations and Sizes of SVC Units Based on Objective Function of Improving the Minimum Voltage 











Magnitude 
Technique SVC location 
Size of SVC Bus 8 Bus 9 Bus 10 
(MVAR) AIS -23.30 3.88 66.19 
MAIS 64.94 50.59 -67.41 





5. CONCLUSIONS 

In this paper, MAIS optimization technique has been used to achieve two different objective 
functions for determining the optimal size of three SVC units installed in the IEEE 14-Bus Reliability Test 
System (RTS). The locations of the SVC units were firstly determined using the VDI approach. When 
targeting one of objective function, the other objective function is slightly having disadvantage. However, 
focusing on each objective will help the Grid System Operators (GSOs) to compare and decide which 
objective need to be considered depending on the current situation and condition. It can be concluded that 
MAIS is superior than the original AIS for solving static VAR compensator planning problem in terms of 
giving lower total system loss and improving minimum voltage magnitude. This due to multiple stages of the 
same processes of cloning, mutation and selection in MAIS that increase the potential of produced 
individuals in initial population to be selected as the final optimal solution. 

For the future development, the MAIS algorithm can be improved by replacing the Gaussian 
Mutation technique with Lognormal Mutation technique or mutated using other optimization techniques such 
as Particle Swarm Optimization (PSO) or Kinetic Gas Molecule optimization (KGMO) to produce better 
offspring from the mutation process. 
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